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Preface 


In the late 1 990s research was being performed at the Kennedy Space Center to develop in situ 
resource utilization technology for Mars. One study topic was the generation and storage of 
liquid oxygen (LOX) obtained from the atmosphere or regolith, but the transfer of this 
commodity was of concern. Mechanical LOX pumps were deemed potentially too heavy and 
unreliable for an autonomous mission to Mars, and alternatives were sought. One option was to 
use the paramagnetic property of LOX, which is significant enough that small quantities of LOX 
can be suspended against earth gravity with a rare earth magnet. With this application in mind, a 
small, internally funded project was initiated at the Kennedy Space Center in 2000 to study the 
use of pulsed magnetic fields to pump LOX. 

Proof-of-concept testing verified the LOX pumping predictions and resulted in a journal 
publication [1]. Numerous small coils were fabricated on cryogenic flow lines and used to 
produce intense, short-duration magnetic fields resulting in dramatic motion of the LOX. In 
addition, effort was expended on modeling the paramagnetic forces in the LOX, which required 
modeling the magnetic field generated by the coils and the coil inductance, allowing the current 
versus time to be predicted and compared to experiment. 

While we were modeling the coil magnetic field, using Mathematica , we came upon a set of 
simple analytical expressions for the magnetic field and its spatial derivatives in Cartesian, 
cylindrical, and spherical coordinates generated by a simple, infinitesimally thin current loop. 

We wrote a short manuscript documenting these expressions, but did not proceed with 
publication. Instead, this manuscript entered the public domain through two different routes. It 
appeared first on the NASA Technical Reports Server (NTRS) at the following link: 

http://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa. gov/200 1 003 8494 200 1 057024.pdf 

and then later on the Open Channel Foundation, an organization that publishes software from 
academic and research institutions. The link for this second site is given below. 

http://www.openchannelfoundation.org/proiect/view docs.php?group id=288 

Through this link, our manuscript was downloaded multiple times and cited seven times as of the 
writing of this NASA Technical Memorandum, but only by title and author list and sometimes 
by web site URL. We, the authors, have decided that there should be a clear, long-term citation 
for this work, yet we consider it too late to seek formal publication. So we have decided to issue 
a NASA Technical Memorandum, with the same title and authors as the original manuscript, 
allowing it to be located by the search engines, and providing a reliable, long-term source and 
citation for this work. The original manuscript, reformatted, but not in any other way altered, 
appears on the subsequent pages of this memorandum. 

We, the authors, would like to express our gratitude to J.M. Griffith, who found a typo in 
Equation (10) of the original manuscript. This memorandum contains the corrected equation. 

[1] Youngquist, Robert C., Immer, Christopher D., Lane, John E., and Simpson, James C., “Dynamics of a Finite 
Liquid Oxygen Column in a Pulsed Magnetic Field,” IEEE Transactions on Magnetics, 39(4), July 4, 2003, 
pp. 2068-2073. 
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1 INTRODUCTION 

Analytic expressions for the magnetic induction (magnetic flux density, B) of a simple planar 
circular current loop have been published in Cartesian and cylindrical coordinates [1,2], and are 
also known implicitly in spherical coordinates [3]. In this paper, we present explicit analytic 
expressions for B and its spatial derivatives in Cartesian, cylindrical, and spherical coordinates 
for a filamentary current loop. These results were obtained with extensive use of Mathematic a™ 
and are exact throughout all space outside of the conductor. The field expressions reduce to the 
well-known limiting cases and satisfy V B = 0 and V x B = 0 outside the conductor. 

These results are general and applicable to any model using filamentary circular current loops. 
Solenoids of arbitrary size may be easily modeled by approximating the total magnetic induction 
as the sum of those for the individual loops. The inclusion of the spatial derivatives expands their 
utility to magnetohydrodynamics where the derivatives are required. 

The equations can be coded into any high-level programming language. It is necessary to 
numerically evaluate complete elliptic integrals of the first and second kind, but this capability is 
now available with most programming packages. 

2 SPHERICAL COORDINATES 

We start with spherical coordinates because this is the system usually used in the standard texts. 
The Cartesian and cylindrical results in Sections 3 and 4 were derived from the spherical 
coordinate results. 

The current loop has radius a , is located in the x-y plane, is_centered at the origin, and carries a 
current I which is positive as shown (Figure 1 .). It is assumed that the cross section of the 
conductor is negligible. 
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The vector potential, A, of the loop is given by [3]: 


A(r,e)^C' , , cos *’'*' 

Ja 2 + r 2 -2arsin#cos^' 


_/*o 


4/a 


4^ >/a 2 + r 2 + 2arsin# _ 


(2-k 2 )K(k 2 )-2E(k 2 ) 


( 1 ) 


where r, and (p are the usual spherical coordinates, and the argument of the elliptic integrals is 

4«rsin 6 


k 2 = 


a 2 + r 2 + 2arsin6 


( 2 ) 


Note that we use A' 2 for the argument of the elliptic integrals. This choice is consistent with the 
convention of Abramowitz and Stegun [4] where m = Jc. 

For a static field with constant current, the B components in spherical coordinates are [3]: 

B, =—^^(sin6A) (3) 

rsmf) dd 

B „=-~ ( r A v> ) (4) 

B= 0. (5) 


Analytic expression for the field components and their derivatives in spherical coordinates are 
given below. For simplicity we use the following substitutions: a 2 = a 2 + r 2 -2ar sin 6, 
p 2 = a 2 + r 2 + 2 ar sin 6, A 2 = 1 - a 2 //? 2 , and C = p o I/n. We note that if desired, further 
simplifications are possible using various substitutions and groupings. 

Field Components: 




a 2 p 


( 6 ) 


B " = T 20 ■ n U r2 + ° 2 cos 26)E(k 2 )- a 2 K(k 2 )] (7) 

2a psm 6 L J 

Spatial Derivatives of the Field Components: 

~T J ~ = _ 7r 4 -6a 2 r 2 cos 20 1 £ (A 2 )+ 

dr 2ra 4 ft 2 U J V ; (8) 

a’ )]*(»>)} 
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dB, 

dO 


•{ [a 4 -la 


2 r 2 + r 4 + 


_ -Ca 2 

~ 4 sin (9a V 3 

cos 2 6 (-3(7 4 + 2 a 2 r 2 - 3 r 4 ) + 
a 2 r 2 cos4#j£(A 2 ) + 

[2a 2 (a 2 + /- 2 )cos 2 6 J K (A 2 ) | 


dB. 


0 _ 


-C 


r {[«‘-3„ 


v + <*v + 


dr 4a 4 /?Vsin# 

2r b + a 2 (3r 2 - a 2 )(a 2 + r 2 )cos 26 + 
3aVcos40]£(A 2 ) + 


[a 2 (-a 4 + a 2 r 2 -2r 4 + 
a 1 (a 2 - 3/* 2 ,)cos 2#)J (A 2 ) J 


dB H _ -C cos# 

~df 


■{ [ 50 ' 


+ 3a 4 r 2 - 3 a 2 r 


2 „4 


4 a 4 /? 3 sin 2 6 
+2r 6 +(-3a 6 + 2a 4 r 2 +9a 2 r 4 )cos20 + 
a 4 r 2 cos 4#] E(k 2 ) — 

[3a 6 + 2 a 4 r 2 + a 2 r 4 + 2 r b + 

a 2 ( 5r 2 - a 2 )(a 2 + r 2 )cos26 + 

(-la 5 r + la 3 r i - 4 ar 5 ,)sin# + 

aV(a 2 -5r 2 Jsin3#]/f(A 2 ;} 


( 9 ) 


( 10 ) 


(ID 


3 CARTESIAN COORDINATES 

The field components and their derivatives in Cartesian coordinates are given below. These are 
easier to use when rotations or translations are needed and obviate the need to transform the basis 
vectors. The following substitutions are used for simplicity: p =x + y ,r =x +y +z , 
a 2 = a 2 + r 2 - 2 ap, pt 2 = a 2 + r 2 + 2ap, A 2 = 1 - a 2 //? 2 , y = x 2 -y 2 , and C = po I/x. Note that p and 
r are non-negative. 

Field Components: 


B x = [(a 2 + r 2 )E(k 2 )- a 2 K(k 2 )] (12) 

B y = J^r[( a2 +r 2 )E(k 2 )-a 2 K(k 2 )]=^ B x ( 1 3) 
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B z = ^-[(a 2 -r 2 )E(k 2 ) + a 2 K(k 2 )] 

Spatial Derivatives of the Field Components: 

it [“ 4 (~ y<3z2 + -y 2> )- 

a 2 (p 4 (5x 2 + y 2 )-2p 2 z 2 (2x 2 + y 2 ) + 3z 4 y)~ 
r 4 (2.x 4 + y(y 2 + z 2 ,))J E (k 2 ) + 

| a 2 (y(a 2 + 2z 2 )- p 2 (3x 2 -2y 2 )} + 
r 2 (2x 4 + y(y 2 + z 2 )^j a 2 A^(& 2 )J 

{[*<*-* 

2 a 6 - 2a 2 (2p 4 - p 2 z 2 + 3z 4 )] E(k 2 ) + 

[r 2 (2r 2 + p 2 )-a 2 (5p 2 -4z 2 )+2a 4 ^a 2 K(k 2 )} 

2z 2 (a 4 -6a 2 p 2 + p 4 )+ z*(a 2 + p 2 )]e(i< 2 )- 
[p 2 -a 2 ) 2 + z 2 (p 2 + a 2 )]a 2 K{k 2 )} 

MK _ dB x 
dx 3y 


(14) 


(15) 


(16) 


(17) 

(18) 




a 2 (p 4 (5y 2 +x 2 )-2p 2 z 2 (2y 2 +x 2 )-3z 4 y}- 
r 4 (2y 4 -y(x 2 +z 2 ))^E(k 2 ) + 

\a 2 (~y(a 2 +2 z 2 )- p 2 (3 y 2 -2 x 2 j) + 
r 2 (2_y 4 - y(x 2 + z 2 a 2 K (k 2 ) J 

(19) 

^ v §0 
1! 

X s- 

(20) 

s1<« 

II 

si* 

(21) 
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BB : _ BB, 
By Bz 


lt = y^¥ a 2 <f ,2 - z '>- la,+r '} E<k2 > + 

a 2 [a 2 -r 2 ]K(k 2 ) 


( 22 ) 


(23) 


4 CYLINDRICAL COORDINATES 

9 9 9 9 

The following substitutions are used for simplicity: a=a+p+z- lap, 
ft 2 = a 2 + p 2 + z 2 + lap, k 2 = \ - a 2 //} 2 , and C = po I/n. 

Field Components: 

B r = + P 2 + z2 )E(k 2 )- a 2 K(k 2 )] 

2 a pp'- J 

c 


B z = " P 2 ~ z " )E(k 2 )+a 2 K(k 2 )] 


(24) 

(25) 


Spatial Derivatives of the Field Components: 


a 4 (3z 2 -8 p 2 ) + a 2 (5 p 4 -4 p 2 z 2 + 

?>z 4 f\E (k 2 )-a 2 [a 4 -3a 2 p 2 + 2p 4 + 

(2a 2 + 7>p 2 )z 2 + z 4 ^K(k 2 )) 

2z 2 (a 4 -6a 2 p 2 + p 4 j\E(k 2 )- 

a 2 [(a 2 - p 2 ) 2 +(a 2 + p 2 )z 2 ]K(k 2 )} 


BB, Cz 


3z 2a^{[( 6<,2f/,J - z2j - 7a,+ 

(p 2 + z 2 ) 2 \E(k 2 )+ a 2 [a 2 - p 2 - z 2 ] K (k 2 ) } 


(26) 


(27) 


(28) 


BB. _ dB lt 
dp Bz 


(29) 
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5 LIMITING CASES 

Several special limiting cases are given for completeness. We have confirmed that our results 
given above do indeed converge to these fonnulas. We also give additional expressions for B x 
and B v near the axis that may prove useful. 

Along the axis of the loop: 


B 


2 


2 (a 2 + z 2 )% 


Near the axis of the loop ( x,y«a ): 

_ 3 a\lxz 
X 4 (a 2 +z 2 /> 

B _ 3a 2 / ti 0 /yz 
' 4 (a 2 + z 2 )% 


Far from the loop ( r»a ): 

2 n r 

d Bo /i 2 vSin# 

b h =-j L (Ina ) — — 
An r 


(30) 

(31) 

(32) 

(33) 

(34) 


6 CONCLUSIONS 

We have presented simple, closed-form algebraic formulas for the magnetic induction and its 
spatial derivatives of a filamentary current loop that are exact everywhere in space outside the 
conductor. Although these formulas are exact, they do require the numerical evaluation of 
elliptic integrals. 

Solenoids with circular cross sections of arbitrary size and configuration can be modeled by 
simply summing the contributions of each individual loop. 

There are, of course, other ways to obtain B for the basic circular current loop. For example, 
series expansions are available [3] or numerical integration via a finite element approach can be 
performed. However, these suffer from limitations such as truncating the series expansions after 
some tolerance is reached or accepting some graininess when using a discrete grid. Our approach 
has neither of these limitations and yields results are that exact up to the limitations of the 
numerical arithmetic and the elliptic integral routines. 

The inclusion of the spatial derivatives allows convective derivatives to be found exactly and 
may prove useful for magnetohydrodynamics applications. 
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